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Abstract. The kinetic features of plasmoid chain formation 
and evolution are investigated by two dimensional Particle- 
in-Cell simulations. Magnetic reconnection is initiated in 
multiple X points by the tearing instability. Plasmoids form 
and grow in size by continuously coalescing. Each chain 
plasmoid exhibits a strong out-of plane core magnetic field 
and an out-of-plane electron current that drives the coalesc- 
ing process. The disappearance of the X points in the co- 
alescence process are due to anti-reconnection, a magnetic 
reconnection where the plasma inflow and outflow are re- 
versed with respect to the original reconnection flow pattern. 
Anti-reconnection is characterized by the Hall magnetic field 
quadrupole signature. Two new kinetic features, not reported 
by previous studies of plasmoid chain evolution, are here 
revealed. First, intense electric fields develop in-plane nor- 
mally to the separatrices and drive the ion dynamics in the 
plasmoids. Second, several bipolar electric field structures 
are localized in proximity of the plasmoid chain. The anal- 
ysis of the electron distribution function and phase space re- 
veals the presence of counter-streaming electron beams, un- 
stable to the two stream instability, and phase space electron 
holes along the reconnection separatrices. 

Keywords. Magnetic Reconnection, Plasmoid Chain, Plas- 
moid Coalescence, Core Magnetic Field, Two-stream Insta- 
bility, Electron Holes 



1 Introduction 

Magnetic reconnection is believed to be the key engine to 
conversion of magnetic energy into heating and accelera- 
tion of space and astrophysical plasmas. Energy is typically 
stored as magnetic energy, then suddenly released released 
in form of violent plasma jets and non-thermal particle ac- 



celeration. Magnetic reconnection is a basic plasma physics 
phenomenon triggering space weather events such as magne- 
tospheric substorms or solar flares. Because of the ubiquity 
of this phenomenon in nature and because it is at the base 
of macroscopic solar physics and magnetospheric events, 
many theoretical, computational and experimental studies 
have been devoted to explain the origin and the dynamics 
of magnetic reconnection. However, one main question has 
still no certain answer. It is unclear why magnetic reconnec- 
tion in nature proceeds on faster time scale than the resistive 
diffusion timescale predicted by the Sweet-Parker model in 
the magnetohydrodynamics (MHD) framework. It has been 
shown the all simulation models including the Hall physics, 
based the decoupling of electrons and ions at ion skin depth 
scale, enable reconnection to occur over faster time period 
{Bim et al. 2001[ l. It has been proposed recently that recon- 
nection occurring in multiple points and secondary plasmoid 
generation can lead to fast reconnection in the MHD con- 



text also via stochastic plasmoid chain formation {Lazarian 



and Vishniac . 1999 ; Kowal et al. . 2009 ; Lapenta', 2008 Bhat- 
(tacharjee et al.^ 2009 ; Uzdensky et al.^ 2010^) . The present 
paper studies the combination of these two mechanisms to 
achieve fast reconnection: a kinetic model of plasma, where 
Hall physics is a built-in feature, is used to study plasmoid 
chain magnetic reconnection. 

The goal of this paper is to study the formation and evolu- 
tion of plasmoid chain structures in the kinetic framework by 
two dimensional Particle-in-Cell simulations. Plasmoids are 
high density structures forming from the outflow plasma ac- 
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celerated in magnetic reconnection {Hones 1977) . A closed 
magnetic loop develops in correspondence of the plasmoids, 
this is why the expression magnetic island is often used in- 
stead of plasmoid. The boundary layer between the plas- 
moid and the inflow plasma is called separatrix. If magnetic 
reconnection occurs in multiple points, the outflow plasmas 
from adjacent reconnection sites form multiple plasmoids, 
organized as beads in a chain. Plasmoid chains are highly 
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dynamic configurations: plasmoids can coalescence forming 



larger plasmoids (Bhattacharjee et al. 1983 1, bounce {Intra 
tor et al. 2009'; Karimabadi et al. 201 1 1 and eject smaller 



plasmoids (haughton et al. 2009| l. Observational studies of 
solar and magnetospheric plasmas confirmed the presence 
of plasmoids (flux ropes in three dimensions) in magneto- 
spheric and solar environments. The plasmoid chains have 
been observed in sun corona i\Bemporad . 2008| and in Earth 
magnetopause and magnetotail ( |Liw et al.\ |2008t \KhotyainF\ 
\sev ef a/.|[20T0l l. 

The formation and evolution of plasmoid chain are a topic 
of many recent theoretical and computational investigations. 
Early simulations studied the occurrence of plasmoid chain 



by MHD simulations {Biskamp 1986 Malara et al. 1992 1. 
Theoretical models and simulations showed how the occur- 
rence of secondary plasmoids affects magnetic reconnec- 
tion and and how fast reconnection in low resistivity plas- 
mas can be achieved in the MHD framework via stochastic 



magnetic reconnection ( Lazarian and Vishniac 1 999 ; Kowal 
eFal] |2009l \Loureiro et al.\ |2005| |2007[ \Lapenta\ |2008 



Bhattacharjee et al. 2009t Cassak et a/.[|2009[ Uzdensky\ 



et al. 2010) . The realization of fast magnetic reconnection 



has also been investigated by Hall MHD simulations [Baal 



\rud et al. 201 l| l. Collisionless Particle-in-Cell simulations 
\Tajima et oT] [7987} \Pritchett\ [20081 P^a al\ |2010b|a 



Tanaka et al. 2010| l studied the dynamics of plasmoid chain. 



The focus of previous Particle-in-Cell studies was primarily 
on the mechanisms that lead electron acceleration to rela- 
tivistic energies. Differently from previous works on colli- 
sionless plasmoid chain, the goal of this paper is to study the 
general properties of magnetic reconnection in a plasmoid 
chain in the kinetic framework, giving a broad picture of the 
plasmoid formation and evolution. Two new results are here 
reported. First, an intense electrostatic activity driving the 
ion dynamics is detected along the reconnection separatrices 
in the plasmoid chain. Second, counter-streaming electron 
beams along the separatrices are unstable to the two-stream 
instability and electron phase space holes appear as a result 
of this instability. 

The paper is organized as follows. First, the simulation 
model is presented. Second, the formation and evolution 
of plasmoid chain, the magnetic and electric field configu- 
rations are analyzed; a study of ion and electron distribu- 
tion functions and phase space are presented to unveil pos- 
sible acceleration mechanisms and the presence of micro- 
instabilities in proximity of the plasmoid chains. Third, the 
simulation results are discussed and compared with those of 
previous works. Finally, the results are summarized and fu- 
ture work is outlined in the conclusion section. 



2 Simulation set-up 

Particle-in-Cell simulations are carried out in a two dimen- 
sional space, the X ~Y plane with Z as ignorable direction. 



with three components for particle velocities and field com- 
ponents. The X, y coordinates correspond to the —xgsm 
(Earth-Sun direction in magnetotail), zgsm (North-South di- 
rection in magnetotail) coordinates in the Geocentric Solar 
Magnetospheric (GSM) system. 

A double current sheet with reversing magnetic field in the 
X direction is initialized as: 



BJB, ^ tanla( ^'^/^ - tanh( " ^^^/^ ^ 



1. 



(1) 



A ' 'A 
A uniform guide field is present along the out-of-plane Z 
direction and equal to Bq/2. The initial density profile is a 
double Harris current sheet with a superimposed background 
plasma with density nt: 



n{y) - 0.2 no f cosh-^ (ILJ^I) +cosh-2 (^LJMI; 

V A A ^ 

where no is the reference density and n;, = 0.2 no; L^ — 
200 di, Ly = 30 di and A = 0.5 di are the simulation box 
lengths in the X and Y directions and the half-width of the 
double current sheet respectively. The lengths are normalized 
to di — c/utpi, the ion inertial length, with c the speed of light 
in vacuum and ujpi = -y/ (47rnoe^/mi). The ratio between the 
ion plasma and cyclotron frequencies is LOpi/fld — 1030.9. 
The ion mass rn^ is taken 256 larger than electron mass me 
and e is the elementary charge. The plasma is initialized 
with a Maxwellian velocity distribution. The electron ther- 
mal velocity is chosen vthe/c — (T'e/(meC^))^/^ = 0.0045. 
Ion temperature is = 5 Tg. The simulation parameters are 
chosen to mimic the plasma characteristics in the solar wind 
current sheets as observed by Wind and Cluster spacecrafts 
{Eri ksson and al. \ |2012| l. In order to speed-up the occur- 
rence of magnetic reconnection, instead of using a perturba- 
tion as in previous studies, the initial system configuration 
is chosen not in equilibrium. The initial current is 20% of 
the current necessary to support consistently the initial mag- 
netic field configuration. As result of this non-equilibrium, 
the plasma is initially accelerated toward the double current 
sheet to establish a current consistent with the initial mag- 
netic field configuration, accelerating the reconnection initi- 
ation. 

The simulation time step is cUpiAt = 0.3. The grid is 
composed of 2560 x 384 cells, resulting in a grid spacing 
Ax = Ay = 0.078 d.^. In total 7 x 10^ computational parti- 
cles are in use. Periodic boundary conditions for particles 
and fields are applied in the X and Y directions. Simu- 
lations are carried out with the parallel implicit Particle-in- 
Cell iPIC3D code (Markidis et a/.||2010| l. These simulations 
would be extremely difficult for be carried out with explicit 
Particle-in-Cell codes because of the size of the simulation 
box and the simulated time period. If compared with typi- 
cal parameters in explicit Particle-in-Cell simulations, large 
time step, grid spacing and ion to electron mass ratio have 
been used. Numerical stability is retained, thanks of the im- 
plicit discretization in time of the Particle-in-Cell govern- 



,(2) 



ing equations {Lapenta et al. 2006 1. In fact, a time step 
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At = 0.3 ujp^^ — 4.8 ujp^ makes the explicit Particle-in-Cell 
method unstable in these parameters regime. Moreover, the 
grid spacing in use is approximately 277 Debye lengths. The 
possibility of using grid spacings that are much larger than 
Debye length still avoiding aliasing instabilities, such as the 
finite-grid instability, is a property of the implicit Particle- 
in-Cell. Moreover, an ion to electron mass ratio equal to 
256 is here used, instead of 25 and 100 in previous stud- 
ies with explicit Particle-in-Cell codes i\Pritchett\ |2008[ \Oka\ 
^TaL\ |2010b|a| \Tanaka et aL\ |20To|. The effect of the im- 
pUcit numerical scheme is to damp artificially the unresolved 



waves, such as the light and Langmuir waves iBrackbill and 



Forslund 1985 Markidis and Lapenta 201 \) . The implicit 
Particle-in-Cell algorithm still correctly describes phenom- 
ena evolving on time scales that are larger than the time step. 
Magnetic reconnection occurs over time period of tens of 
f2~^, approximately ten thousand times the simulation time 
step, and the implicit numerical scheme does not alter its evo- 
lution. 



3 Plasmoid chain formation and evolution 

A tearing instability rapidly occurs out of particle noise and 
generates several plasmoids. The reconnected flux for the 
single plasmoid ijj is calculated as the difference of the out- 
of-plane component of the vector potential at X and O 
points along the line y — Ly = 22.5 di (center of the top 
current sheet). Figure 1 panel a shows a stack plot of the 
out-of -plane component of the vector potential A^, normal- 
ized using the background ion skin depth and the asymptotic 
magnetic field Bq. An A^ minimum identifies an X point 
(green and red colors in Figure 1 panel a, while a maxi- 
mum (black color) indicates an O point. The plot illustrates 
multiple plasmoids collapse, forming larger plasmoids. The 
reconnection rate is used to characterize the speed at which 
magnetic reconnection occurs. A separate analysis of the re- 
connection rate at different reconnection sites is required, be- 
cause magnetic reconnection develops in a large number of 
points. In addition, the reconnection sites move in time, it 
is therefore necessary to identify the reconnection site and 
follow them during the simulation. The reconnection rate 
is calculated for each reconnection site, indicated with solid 
lines in Figure 1 panel a. The reconnection rate can be com- 
puted in two ways: (i) the time derivative of the reconnected 
flux tp (solid lines in Figure 1 panel b) (ii) the reconnection 
electric field in the upstream region, normalized to the 
Alfven velocity Va and Bq (dashed lines in Figure 1 panel 
b). The two measures of reconnection rates are in agreement 
in each site. The reconnection rate evolves very similarly for 
aU the X points: it grows very rapidly until time i^cit = 1-5, 
when it reaches the 0.3 peak value, then decreases reaching 
the range 0.01 — 0.08 for different X points. 

As result of magnetic reconnection, several high density 
structures are formed by the outflow plasmas of adjacent 



X points (Figure 2, time ildt = 3.5). The electron plasmoid 
density is approximately four times the background density. 
At time fldt — 3.5, plasmoids are large 5 di and progres- 
sively grow by merging and forming larger plasmoids. Low 
density layers, called cavities iPritchett 2008| l, form along 
the reconnection separatrices. During the simulation, the ex- 
pulsion and formation of smaller plasmoids are not observed. 
Because secondary islands are not produced during magnetic 
reconnection, the plasmoids have approximately the same 
spatial size and very large plasmoids, such as monster plas- 



moids ( Uzdensky et al. 2010 1, do not appear. The hierarchi- 
cal self-similar structure of the plasmoid-dominated current 
sheet is not detected in the Particle-in-Cell simulations ( |t/z- 
\densky et a/. [ [20 10). 

All the plasmoids develop a relatively strong out-of -plane 
core magnetic field, as shown in the B^ contour plot in Fig- 
ure 3. The core magnetic field is unipolar and in the same 
direction for the plasmoids of both current sheets. The inten- 
sity of Bz magnetic field without the guide field is approxi- 
mately 0.5 Bq, that is the value of the guide field in the pro- 
posed simulation. As noted in previous works ^Karimabadi 
et al. 1999| l, the core magnetic field is not uniform over the 



plasmoids, but it presents a complex pattern. The S-shaped 
Bz regions, that were found in hybrid simulations of multi- 
ple X points reconnection, are present also in the B^ contour 
plot in Figure 4. The coalescence process drives the anti- 
reconnection I^Oka et al.\\2QlQhl\Tanaka et al.\\2QlQ) , where 
magnetic reconnection inflow and outflow are reverted. The 
inflow is directed in the X direction, while the outflow direc- 
tion in anti-reconnection is along the Y direction (white ar- 
rows in Figure 4). The signature of anti-reconnection is vis- 
ible in the typical quadrupolar structure of the Hall magnetic 
field in the white dashed squares in Figure 4. The quadrupole 
Hall magnetic field is rotated 90° with respect to the typical 
reconnection Hall magnetic field and generates the S-shaped 
regions, observed previously in multiple X points reconnec- 
tion hybrid simulations ^Karimabadi et al. 1999| l. 

The out-of-plane current component of the current drives 
the coalescence process. In fact, plasmoids act as wires car- 
rying current in the same direction attracting each other, co- 
alescing one plasmoid into another one. The total current 
is supported only by the electrons: the intensity of electron 
current is seven times higher than the ion current, developing 
in the opposite direction (Figure 5). This result is in agree- 
ment with the three-dimensional kinetic simulations of flux 
ropes, where the axial current is supported only by electrons 
\Markidis et a/.[[20lTal l. 

The most intense electric fields are in-plane and local- 
ized along the separatrices or in the plasmoids, as shown 
in Figure 6. The maximum value of the electric field is 
~ 2 VaBq/c, that is approximately 10 times the peak recon- 
nection rate. 

The electric field configuration in a plasmoid region is 
shown in detail in an enlargement of an area occupied by 
a plasmoid at time Vlcit = 18.62 in Figure 7. A strong inward 
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Fig. 1. The stack plot of the out-of-plane vector potential component Az along the line y = 22.5 di is shown in panel a. The reconnection 
rate in panel b) is calculated as derivative in time of the reconnected flux (solid line) and as Ez (dashed lines), for different X points 
(colored lines in panel a. 



electric field is present in its center. Along the separatrices, 
the electric field is normal to the separatrices (black line in 
panel b): it is directed inward inside the separatrix, while 
it is outward just outside the separatrix. Thus, ions in the 
plasmoid are accelerated inward, while electrons are accel- 
erated outward and reflected back in the plasmoid. The ion 
bulk flow in a plasmoid is presented in Figure 7 panel b. In 
guide field reconnection, the reconnection jets are not com- 
pletely directed along the X direction, but their directions 



are tilted by the guide field ( Goldman et al. 201 1 1. Thus, the 



reconnection jets from different X points do not run head- 
to-head, but stream one over the other. The strong ion bulk 
flow is localized inside the plasmoid. An analysis of the ion 
phase space Vx and Vy vs x in a 0.234 di thick region along 
y = 22.5 di (blue dashed line in panel b) shows that ions are 
accelerated towards the plasmoid center, and in the — Y and 
Y directions for small and large X values respectively. This 
combined motion generates a rotation of the ion bulk flow in 
the counterclockwise direction. Outside the separatrix line, 
a weak ion flow rotates in the opposite direction. This ion 
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Fig. 2. Contour plots of the electron density, normalized to the background density nt, at different times. The plot shows the formation and 
the progressive growth of relatively high density regions, the plasmoids, by coalescence. 



(Bz-Bg)/B0 




Fig. 3. Contour plots of the out-of-plane magnetic field without the guide field, B^ — Bg, normalized to Bo at different times. Each magnetic 
island is characterized by an out-of-plane core magnetic field, whose value is approximately equal to the guide field Bg = 0.5-Bo. 
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Fig. 4. Contour plots of the out-of-plane magnetic field without the guide field, — Bg, normalized to Bq at time fldt = 26.77. The 
core magnetic field shows a complex pattern. The S-shaped regions are visible in the coalescing area. The quadrupole structure of the Hall 
magnetic field in the coalescence-driven anti-reconnection is enclosed in white dashed rectangles. The white arrows indicate the outflow 
direction in anti-reconnection. 
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Fig. 5. Contour plot of the out-of plane electron and ion currents, J^e and Jzi, normalized to enoc, at time Qcit — 18.62. The total current 
is supported only by the electrons while the the ions generate a current in the opposite direction. 



bulk velocity system is consistent with the ion cuiTent pat- 
tern, presented in Figure 5 panel b. 

In magnetic reconnection studies, the bipolar electric field 
signatures are regarded as flag of electron micro-instabilities 
{Lapenta et al. 201 1 1. An investigation of the parallel elec- 
tric field i?// = E-B/|_B| in Figure 8 reveals that there 
are several bipolar parallel electiic field structures inside 
the plasmoids and along the separatrices. These signatures 
have maximum and minimum values ^Q.ZBqVa/ c, approx- 



imately equal to peak reconnection electiic field. 

An enlargement of the same region presented in Figure 7 
shows that bipolar electric field structures are aligned along 
the magnetic field lines, suggesting the presence of a beam 
instability. The electron distribution function is analyzed to 
verify this hypothesis. The electron distribution functions, 
fe{vx,Vy), in four different locations are presented in Fig- 
ure 9 panel a. The four different regions, black boxes in the 
central subpanel, have been chosen along the separatrix (red 
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line in the E j j contour plot in the central subpanel), covering 
a 0.234 di x 0.234 di wide area. A comparative analysis of 
the four electron distribution functions shows the presence of 
counter-streaming electron beams along the magnetic field. 
In subpanel I, the electron distribution is taken in proximity 
of an X point and is approximately Maxwellian. In subpanel 
II, the distribution function shows an electron core popula- 
tion at positive Vx and Vy, a beam moving from left to right 
along the magnetic field direction (red line in the subpanels). 
An electron beam at negative and Vy is present also. Sub- 
panels III and IV confirm the presence of two beams, one at 
positive Vx and Vy and one at negative Vx and Vy, moving in 
opposite directions. 

Figure 9 panel b shows the electron phase space along a 
separatrix, x// vs u// (curvilinear coordinate along the mag- 
netic field line vs the parallel velocity). For high and low a;// 
values, two electron beams are visible as two electron pop- 
ulations moving at opposite speeds. The interaction of two 
beams results in the formation of phase space electron holes, 
in white dashed lines in Figure 9 panel b. There is clear 
indication that electron holes are generated by two-stream 
instabilities of the two electron beams along the separatrices. 



4 Discussion 

Magnetic reconnection in a plasmoid chain has been stud- 
ied by kinetic simulations. Magnetic reconnection is initially 
triggered by the tearing instability in multiple X points. Sev- 
eral plasmoids originate from the outflow plasma of magnetic 
reconnection. The reconnection rate calculated for different 
points reaches a maximum in a very short period fici^ ~ 1-5. 
After a very fast transient, the plasmoid chain dynamics is 
dominated by plasmoid coalescence. In this second phase, 
the reconnection rate decreases rapidly reaching an asymp- 
totic values range 0.01 — 0.08 for different X points. 



Several high density regions, the plasmoids develop as a 
result of the initial tearing instability. Plasmoids grow in size 
by merging. At a given time, the plasmoid are approximately 
the same size. Secondary plasmoid formation and ejection 
of smaller plasmoids are not observed. The self-similar plas- 
moid hierarchy in the plasmoid chain (Uzdensky etal. 2010| l 
is not present in the reported simulation. This is in contrast 
with previous plasmoid chain Particle-in-Cell simulation re- 
sults ( Oka et al. 2010b|a l), where secondary plasmoids were 
indeed detected. The difference between previous Particle- 
in-Cell simulations with the here presented results, might 
be explained by the use of different simulation box lengths, 
plasma temperatures, and ion to electron mass ratios. It is 
important to note that collisionless Particle-in-Cell simula- 
tions with realistic mass ratio do not report the formation of 
secondary plasmoids ^Lapenta et al. 2010[ Y^apenta et al. 



201 1 ). Further studies are necessary to identify which sim- 
ulation parameters determine the occurrence of secondary 
plasmoids in Particle-in-Cell simulations. 

An out-of-plane core magnetic field develops in corre- 
spondence of the plasmoids. The strong core magnetic field 
has been observed several times in satellite observations. 
ISSE 3 and GEOTAIL missions reported the presence of 



cross-tail magnetic field in the Earth's magnetotail iSlavin 



eToI] [19891 1. M oreover, MHD ^Hesse et al.\ \i996} and hy- 
brid simulations ( Karimabadi et al.\ 1999| l confirmed the for- 
mation of a plasmoid core field and proposed a formation 
mechanisms in terms of Hall-generated currents. Moreover, 



hybrid simulations {Karimabadi et al. 19991 revealed that 



the core magnetic field has a complex structure originating 
from the complex flow pattern in the plasmoid. In multi- 
ple X points reconnection hybrid simulation, S-shaped 
regions were found. The fully kinetic simulation in this pa- 
per confirms the formation of unipolar core magnetic field 
and complex pattern. Some of these structures, such as 
S-shaped regions, can be explained as a result of the Hall 
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magnetic field in anti-reconnection. 

The most intense electric fields develop in-plane, reaching 
a peak value that is approximately ten times the maximum 
reconnection electric field. They are localized along the sep- 
aratrices: the electric fields point inward/(resp. outward) in- 
side/(resp. outside) the plasmoid. Intense normal electric 
fields have been observed in magnetotail and an ion accel- 
eration mechanism due by repetitive ion reflections between 



the separatrices have been proposed {Wygant et al. 2005 i. 



The kinetic simulations in this paper shows an electric field 
configuration, that can support such acceleration mechanism 
{Wygant et al. 2005) and is agreement with the results of 
Particle- in-Cell studies ^arkidis et al. \ |201 Vo\c\ . The ki- 
netic study of plasmoid chain reveals that the separatrices 
are regions populated by counter-streaming electron beams, 
eventually driving the two-stream instability. An analysis of 
the electron distribution functions and phase space along the 
separatrix shows the presence of electron beams along the 
separatrices and the formation of electron phase space holes. 
Bipolar electric field structures were previously observed in 
{Drake et al. 20051. In this study, it is proven that bipo- 



lar electric field structures are associated with electron phase 
space holes, produced by the two-stream instability of elec- 
tron beams along the separatrices. The two-stream instability 
could generate the electron holes observed in plasmoid em- 



bedded in the magnetotail current sheet {Khotyaintsev et al. 
|20T0l l. 



5 Conclusions 

The general properties of magnetic reconnection in a plas- 
moid chain in the kinetic framework have been studied by an- 
alyzing two dimensional Particle-in-Cell simulation results. 
A plasmoid chain is rapidly formed by the tearing instability 
and plasmoids coalescence progressively in larger plasmoids. 
Secondary islands are not formed and plasmoids are approx- 
imately the same size with no hierarchical self-similar struc- 
ture. This is probably a result of the chosen simulation pa- 
rameters. Future work will investigate which quantities reg- 
ulate the secondary plasmoid generation in Particle-in-Cell 
simulations. Each plasmoid is characterized by a core out-of- 
plane magnetic field and an out-of -plane electron current that 
drives the coalescence. The quadrupolar structure of the Hall 
magnetic field is the signature of anti-reconnection. Bipolar 
electric field structures are present along the plasmoid sep- 
aratrices. The analysis of the electron distribution functions 
and phase space show counter-streaming electron beams, un- 
stable to the two-stream instability, and electron phase space 
holes along the reconnection separatrices. 

The plasmoid chain formation and evolution have been 
presented in two dimensional simulations. The reduced ge- 
ometry allows to simulate the plasmoid chains with realis- 
tic parameters in relative large systems (hundreds di), but 
the physical phenomena developing in the third direction 



are neglected. In real three dimensional systems, the mag- 
netic islands correspond to extended flux ropes, interact- 
ing in a variety of complex ways non possible in two di- 
mensional geometries. The comparison of two and three 
dimensional MHD simulation results show significant dif- 
ferences in the evolution of turbulent reconnection ^Kulpa- 
\Dybel et al.]pMO \^o wal et al.\^W9\ . In the kinetic frame- 
work, micro-instabilities in the third dimensions, such as the 
lower hybrid instabilities ^Ricci et ai\\2005\ , are not included 
in the model. In fact, recent three-dimensional simulations 



of magnetic reconnection with guide field ( Daughton et al. 



2010 2011 1 show that oblique tearing modes can generate 



turbulence along the reconnection separatrices. Future three 
dimensional studies are necessary to evaluate the effect of 
flux rope interaction and waves propagating in the third di- 
rection on the plasmoid chain evolution. 
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Fig. 7. Enlargement of a simulation region, x — 84 — 93, j/ = 15 — 30 enclosing a plasmoid at time ridt — 18.62. A quiver plot of the electric 
field and ion bulk flow are shown in panels a and b. The length of the arrows indicates the intensity. The red arrows shows the direction of 
the ion bulk flow. The phase-space plots, Vx and Vy vs x al y = 22.5 (blue dashed line in panel b) is shown panel c. 
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Fig. 8. Contour plot of the parallel electric field E j j = E- B/|_B|, normalized to the reconnection electric field Bi^VaIc with superimposed 
magnetic field lines in black at time f2cii = 18.62 and enlargement of of a simulation region, a; = 84 — 93, j/ = 15 — 30 enclosing a plasmoid. 
Bipolar parallel electric field signatures are localized in proximity of the magnetic separatrices. 
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Fig. 9. Electron distribution functions, s// vs u// (curvilinear coordinate along the red line in panel a vs the parallel velocity), at four different 
locations (black boxes) at time Qcit — 18.62 in panel a. The red line in the distribution functions shows the direction of the magnetic field. 
Electron phase space along the red line of panel b with phase space holes, encircled in white dashed lines. The electron distribution functions 
show the presence of counter-streaming electron beams, while the phase space plots reveal the presence of electron phase pace holes along 
the separatrices. 



